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Abstract. We discuss recent improvement in the treatment of gluon radiation in top 
production and decay in e + e~ processes according to the HERWIG event generator 
and show studies on the top mass reconstruction at the future Linear Collider. 
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For the sake of performing precision measurements of top quark properties at 
the future Linear Collider, trustworthy Monte Carlo simulations of multiparton 
radiation in top production and and decay will be essential. According to the 
Ph! standard algorithm of the HERWIG event generator [1] , we shall refer to hereinafter, 

multiple radiation is treated in the soft or collinear approximation and no emission 
is permitted in the so-called 'dead zones', which correspond to hard and large-angle 
parton radiation. The HERWIG algorithm can be improved by applying matrix- 
element corrections: the dead zone is populated by the use of the exact first-order 
matrix element ('hard correction') and the 0(cts) result is used in the already-filled 
region any time an emission is the 'hardest so far' ('soft correction') [2]. 

One of the new features of HERWIG 6 [3] consists of the implementation of 
matrix-element corrections to top decays, which have been shown to have a rel- 
evant effect on jet observables at the threshold for top pair production [4]. As 
pointed out in [5], matrix-element corrections to top production in e + e~ annihi- 
lation, implemented following [6], still needed improvement, since mass effects are 
not systematically included in the dead zone boundary and in the soft correction. 

In Fig. 1 we plot the total and HERWIG phase space for the process e + e~ — > 
q(Pi)q(P2)9(P3)i considering massless quarks and top quarks at ^Js = 500 GeV, in 
terms of the energy fractions X\ = 2p\-q/q 2 and x 2 = 2p 2 -q/q 2 , with q = Pi+P2+P3- 
We see that once we account for mass effects, the dead zone includes both a large- 
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and a small-angle region of the physical phase space, the latter corresponding to the 
neighbourhood of the x\ — x<i — 1 point, which, on the contrary, would be entirely 
inside the HERWIG region if we neglected mf/s terms. Since the soft singularity is 
not completely inside the HERWIG region, the total emission into the dead zone, 
naively calculated, would be infinite. As we did for the top-decay case [4], we 
avoid the soft singularity by setting a cutoff E min on the energy of gluons which are 
radiated in the dead zone and check that phenomenological observables are weakly 
dependent on the value of E min . We choose E min = 2 GeV as the cutoff default 
value. 

We consider e + e~ — > ti processes at y/s = 500 GeV and 1000 GeV and cluster 
final-state partons into three jets by the use of the Durham algorithm [7], assuming 
that both IV's decay leptonically. We set the cuts E T > 10 GeV and AR > 0.7 on 
transverse energy and invariant opening angle of clustered jets. In Fig. 2 we plot 
the distributions of ys, the threshold value of the Durham variable for all events 
to be three-jet-like, according to HERWIG 6.2, the latest public version, and 6.3, 
the new version in progress which will fully include mass effects in matrix-element 
corrections to top production. We investigate the options to fill either the small- 
and large-angle dead zone or only the large-angle region. The impact of the full 
implementation of mf / s effects is a suppression of emission, which is more visible at 
y/s = 1000 GeV, as the radiation in the top-production stage gets more important. 
Filling the small-angle region as well results in more events at intermediate values 
of 1/3. We checked that once the centre-of-mass energy is increased so that terms 
ml/ s are negligible, the 6.3 results reproduce the 6.2 ones. 

Matrix-element corrections to W — > qq' decays in the top decay, not yet included 
in [3], turn out to be a straightforward extension of the corrections to Z — > qq 
processes in the massless approximation m q ^ <C ttlw,z- We found little impact on 
generic jet observables at the Linear Collider, even at the top threshold, where the 
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FIGURE 2. y 3 distributions at V* = 500 GeV (a) and ^ = 1000 GcV (b) according to 
HERWIG 6.2 (dashed line) and 6.3, once we populate either the small- and large-angle (solid) or 
only the large-angle region (dotted) of the dead zone. 



radiation in the production phase is negligible. In fact, the already-existing jets 
associated with the b quarks from the top decay and with the W-decay products, 
even in the soft or collinear approximation, make the detection of hard and large- 
angle gluon radiation in the W decay pretty difficult. Systematic and detailed 
analyses to find out possible variables which might be sensitive to matrix-element 
corrections to W decays are in progress. 

We finally wish to report on studies on the top mass reconstruction in the dilepton 
channel. We consider the Mepton invariant mass and the b energy E^, where 
the b quark is considered together with its gluon radiation, as possible variables 
which may allow a fit of the top mass. Being a Lorentz-invariant observable, the 
mu distribution is independent of the centre-of-mass energy and, as already pointed 
out in [8] for the purpose of hadron collisions, of the hard-scattering process as well. 
In fact, within the statistical Monte Carlo fluctuations, we find the same results for 
different values of \fs. On the contrary, E b is not Lorentz-invariant, hence it will be 
sensitive to the boost from the top rest frame, where the top decay is performed, to 
the laboratory frame. At the threshold for ti production, the dependence of E^ on 
the top mass will be emphasized, as the ti pair is produced almost at rest. Fig. 3 
shows that the m u distribution is shifted towards larger values as the top mass 
is increased. Moreover, the half-maximum width of the Et> distribution shows 
a strong dependence on the top mass at \/s = 370 GeV, since the distribution 
gets narrower as the top mass approaches the threshold value y/s/2. If we try to 
parametrize the relation of the average value (mj<) and in terms of m t , we find 
that the best fits are a straight line for (mu) and a parabola for cr&: 

(m u ) = 0.756 m t - 37.761 GeV ; (1) 
a b = -0.081 m\ + 26.137 m t - 2048.968 GeV . (2) 

Inverting Eqns. (1) and (2) to extract m t , we conclude that if A{rribe) and A<Jb are 
the uncertainties on measurements of the invariant mass and of the half-maximum 




FIGURE 3. (a): Distribution of the invariant mass for m t — 171 GcV (dotted line) and 
m t = 179 GeV (solid), (b): b energy y/s = 370 GeV for m t = 179 (solid), 175 (dashed) and 171 
(dotted) GeV. 



width, they will result in an error Am t rs 1.32 A(m M ) and Am t w 0.35 — 0.65 Acr b , 
where the latter uncertainty refers to the range 171 GeV <m t < 179 GeV. The b 
energy looks therefore to be a quite promising observable with which to extract m t , 
although we would need to know the experimental accuracy on the masurement of 
o"5 in order to estimate the foreseen uncertainty on m t . Furthermore, we expect 
that Eb and will be sensitive to the beam energy smearing, which has not been 
accounted for in the plots of Fig. 3. The implementation of beamsstrahlung in 
HERWIG, via an interface with the CIRCE program [9], is under way. 

In summary, we discussed recent progresses in the implementation of matrix- 
element corrections to the HERWIG simulation of top production and decay at 
the Linear Collider and showed studies on the top mass reconstruction, for which 
purpose the b energy is expected to be an interesting variable for e + e~ collisions 
slightly above the top threshold. 
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